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Abstract 

In this paper, we prove that the superhorizon conservation of the curvature perturbation ( in 

single-field inflation holds as an operator statement. This implies that all (^-correlators are time 
i— i independent at all orders in the loop expansion. Our result follows directly from locality and 

^vq diffeomorphism invariance of the underlying theory. We also explore the relationship between 

the conservation of the single-field consistency relation and the renormalization of composite 

operators. 
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1 Introduction 



Inflationary perturbations are likely to admit a weakly coupled description. In particular, the 
observed near-Gaussianity of the primordial fluctuations suggests that a treatment in terms of 
free fields with computably small corrections is applicable. From this point of view, higher-order 
corrections in perturbation theory (i.e. loop corrections) would appear to be unimportant. On the 
other hand, there are examples in field theory in which loop corrections do affect the qualitative 
behavior of weakly coupled systems. For instance, sometimes the tree level contribution to a given 
quantity vanishes for special reasons (e.g. gauge invariance forbids the decay of the Higgs boson 
into photons at tree level). In these cases, loop corrections are the dominant effect (e.g. the Higgs 
decays into photons at one- loop). Moreover, sometimes the coefficients of loop corrections are 
(naively) infinite (e.g. the electron self-energy is divergent). These cases are more subtle because 
they require us to understand the meaning of the divergences in order to achieve physical results 
(e.g. removing ultraviolet divergences may imply renormalization group flow). 

Time-dependent loop corrections to the primordial curvature perturbation £ would fall into 
both categories. First of all, at tree level £ vanishes outside the horizon [1-5] and therefore 
any loop corrections that generate £ ^ would be the leading effect. Second of all, the putative 
corrections are expected to scale as log a(t) [6] and therefore diverge as we take the scale factor a(t) 
to infinity. Such infrared divergences would have to be understood before reliable predictions 
could be made. In order to sharpen the understanding of inflationary perturbation theory, one 
would therefore like to develop results that hold beyond the tree approximation [6-8] . One of the 
most reliable principles for achieving such results is symmetry. 

The universe contains a lot of symmetry. On large scales and/or early times, the spacetime is 
invariant under spatial translations and rotations. This fixes the metric to be of the Friedmann- 
Robertson- Walker (FRW) form 

ds 2 = -dt 2 + a 2 (t)dx 2 , (1.1) 

where dx 2 represents a maximally symmetric three-space (which we will take to be flat space). 
Small fluctuations around the homogeneous background are close to scale-invariant, suggesting 
additional symmetry in the action for the fluctuations (e.g. a global time-translation symmetry). 
The time evolution of the FRW spacetime (1.1) can be thought of as a spontaneous breaking of 
an even larger symmetry group. This symmetry breaking can be characterized by introducing 
the Goldstone mode tt as a perturbation along the broken symmetry, i.e. a local shift in time 
t + tt(x, t). Just as in the classic example of the chiral Lagrangian for pions, the effective action 
for tt is highly constrained by the non- linearly realized symmetry [9]. This approach has been 
particularly fruitful for describing inflation [10] where the time dependence of the couplings for 
the Goldstone mode are constrained by additional global symmetries. For single-clock inflation, 
the Goldstone mode is directly related to the adiabatic fluctuations that are observed in the 
cosmic microwave background (CMB). Transforming to comoving gauge, the field ir is eaten by 
the metric. Adiabatic fluctuations are now represented by the curvature perturbation £, defined 
as the isotropic scalar perturbation to the three-metric, 

g tj = a 2 (t)e 2 ^'% . (1.2) 
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At late times, £ non-linearly realizes conformal symmetries on the spatial slice [2, 11, 12]. For 
instance, under dilatations, x i— > e A a:, the curvature perturbation transforms by a shift, £ i— )■ £+A. 

At the classical level, the symmetries of £ have been used to derive several important the- 
orems about single-field inflation. For example, Maldacena's consistency relation [13] uses the 
dilatation symmetry to show that a long-wavelength mode is unobservable and therefore can't 
induce observable correlations with short-wavelength modes. Similarly, Weinberg's proof [2] of 
the conservation of £ on superhorizon scales uses the dilatation symmetry as a method for finding 
solutions to the classical equations of motion. One might expect that both statements could be 
promoted to operator statements in a quantum mechanical theory (where the symmetry can be 
expressed as a Ward identity [14]). However, by now there are sufficiently many counterexamples 
to both results that it is clear that neither statement should follow from symmetry alone. 

In this paper, we will prove that £ is indeed conserved at all-loop order in single-clock inflation. 
Using nothing more than locality (which forces commutators to vanish outside of the light-cone), 
we first show that the operator £ must satisfy an equation of the form 

t = f[C), (1.3) 

where /[£] is a functional of £ and its spatial derivatives. We then use symmetry [2, 11] to 
constrain the operators appearing in f[Q. Non-derivative operators are forbidden by the dilata- 
tion symmetry. The remaining operators can be organized according to their scaling behavior as 
a — > oo. (To achieve this, we have to define renormalized composite operators [15, 16].) We will 
find that the leading operators on the right-hand side of eq. (1.3) vanish as a~ 2 . This establishes 
that, in the limit a — > oo (or on superhorizon scales), all (^-correlators are time independent at 
all orders in the loop expansion. 

The outline of the paper is as follows: In Section 2, we review the symmetries of adiabatic 
fluctuations in general FRW cosmologies. We show that the curvature perturbation £ non-linearly 
realizes conformal symmetries. We use these symmetries, in Section 3, to provide an all-orders 
proof for the conservation of £ on superhorizon scales. An essential part of the proof is defining 
a renormalization procedure for composite operators in inflationary spacetimes. We relegate 
a technical discussion of this subtle issue to Appendix A. In Section 4, we comment on the 
relationship between our proof for the conservation of £ and Maldacena's consistency relation. 
We state our conclusions in Section 5. 

2 Symmetries of Adiabatic Fluctuations 
2.1 Non-Linearly Realized Symmetries 

Consider an FRW background with a set of matter fields tp m (t). The time dependence of the 
background spontaneously breaks time diffeomorphisms. Just as in particle physics, we can define 
a Goldstone mode n as a perturbation of the fields along the broken symmetry, i.e. a local shift 
in time. This induces adiabatic fluctuations 

8ip m (x, t) =ij) m [t + tt(x, t)) - $ m {t) . (2.1) 
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An effective theory for the Goldstone mode tt has been constructed in [9, 10] (for related work 
see [17-27]). It is clear that in the case of purely adiabatic fluctuations, the perturbations in the 
matter sector can be gauged away by performing a time diffeomorphism 

t\-+t- ir(x,t) . (2.2) 

The fluctuations are then in the metric only. These metric fluctuations are described most 
conveniently in comoving gauge (also called £-gauge), defined as 

Sif) m = and gi j{x,t) = a 2 (t)e 2 ^'% , (2.3) 

where £ is the curvature perturbation [28, 29]. Perturbations in goo an d gto are related to Q 
through the Einstein equations [13]. For simplicity, we will drop tensor fluctuations throughout, 
but re-introducing them doesn't affect our conclusions. In this gauge, the adiabatic mode is 
characterized by (,{x, t) directly. 

From the form of (2.3), we see that the adiabatic mode is invariant under the following large 
gauge transformations 1 [11] : 

dilatation : x 1— > x = xe X , C( x ) ^ ((&) + ^ > (2-4) 
SCTs : x ^ x = x + 2(6 • x)x - x 2 b , ((x) i-> ((x) + 2b • x , (2.5) 

where SCT stands for special conformal transformation. Notice that £ transforms non-linearly: 
dilatations shift the value of £, while SCTs shift its spatial gradient. Both of these symmetries 
are part of the group of diffeomorphisms under which the theory is invariant. What makes the 
transformations in (2.4) and (2.5) special is the fact that they preserve £-gauge, but are not 
removed by gauge fixing. After gauge fixing, the large gauge transformations therefore remain 
a symmetry of the action. As for any global symmetry, this implies the presence of conserved 
currents: one for the dilatation, J%, and three for the special conformal transformations, Jg C uy 
In the following, we will drop the subscripts whenever an expression applies to both types of 
currents and keep it only when a distinction needs to be made. Current conservation, d^J^ = 0, 
implies the following Ward identity [14] for correlation functions [30] 

id^(j»(x,t)(( yi ,U)---((y n ,U)) = 

n 

= TT,6(t - QS(x - Vt) (Civut*) • • • 6^, **)••• C(Vn>**)> . ( 2 - 6 ) 
i=l 

where <5£ denotes infinitesimal variations of £ under the large gauge transformations 

5 d ( = -l-x-d x ( , (2.7) 
J«C = -2x i - 2x\x ■ d x C) + x 2 d\ ■ (2.8) 

Here, we have introduced an index i to distinguish the three SCTs associated with the three 
components of the vector b. Finally, it is also convenient to define a conserved charge associated 
with each symmetry 

Q= [ d 3 x J° . (2.9) 



x By large gauge transformations we mean gauge transformations that do not vanish at infinity. 
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Formally, this satisfies Q = 0. However, when the symmetry is spontaneously broken, IR diver- 
gences make the value of Q ill-defined. On the other hand, Q remains well-defined in commutators 
with local operators, such as [Q, £], and inside correlation functions. In fact, by integrating the 
Ward identity (2.6) for n = 1, we see that 

i[Q, (] = $(■ (2-10) 
2.2 Symmetries and the Conservation of Zeta 

The presence of the dilatation symmetry has played a crucial role in previous work on the con- 
stancy of C outside the horizon. At a technical level, this connection was implemented most 
directly by Weinberg [2], who used the existence of the large gauge transformation to find two 
physical solutions to the classical equations of motion in any FRW background: one solution is a 
constant and the other decays as a~ 3 . These two solutions correspond to the growing and decay- 
ing contributions of the adiabatic mode. If we assume that only the adiabatic mode is present, 
then we have found all the possible solutions and therefore ( is conserved classically. 

Tree-level. — Although Weinberg used Newtonian gauge, his result is easily reproduced from the 
dilatation symmetry in ("-gauge. For our purposes, it will be useful to state Weinberg's proof in 
a quantum mechanical language using the Ward identity (2.10). Taking the expectation value, 
we find 

([Q d ,(k}) = i(2nf5(k) . (2.11) 
Since Qd = 0, the time derivative of this expression is 

([Qd,(k}) = 0. (2.12) 

In order to satisfy (2.11), we require a non-zero solution for (fc_>-o> while (2.12) implies that this 
solution is time independent. We have therefore found that a non-zero constant is a solution for 
Cfc— • Moreover, locality requires that (see §3.1) 

[Ut), C fe '(*)] « a~ 3 (t) (2ir) 3 5(k + k') . (2.13) 

This implies the existence of a second solution scaling as a~ 3 . Since there are only two solutions 
to the classical equations of motion, we have found that ( is classically conserved. 2 In this paper, 
we will extend Weinberg's proof to the quantum level. 

One-loop. — It is well-known that massless scalar fields can receive time evolution outside the 
horizon from quantum corrections. Essentially, this arises because radiative corrections induce 
a mass for any unprotected scalars, which then sources superhorizon evolution [8]. Two-point 
functions are found to evolve as loga(i). It is therefore natural to ask what happens to the 
conservation of ( at loop level. This question was first raised by Weinberg in [6, 7]. Subsequently, 
a calculation by Kahya, Onemli and Woodard [31], indeed, suggested that loops would induce a 
time dependence of (. This conclusion was challenged by Pimentel, Senatore and Zaldarriaga [32]. 
In an impressively complex calculation, these authors showed that although individual one-loop 

technically speaking, we have not shown that these solutions can be extended to finite momentum k. However, 
using the Ward identity (2.6) it is straightforward to prove that this is the case (see Appendix A of [30]). 
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diagrams do induce a time dependence, the effect precisely cancels when all diagrams are summed. 
Not surprisingly, symmetry played an important role in understanding this cancellation. At 
various stages in their calculation Pimentel, Senatore and Zaldarriaga, directly or indirectly, 
employed the dilatation symmetry. They showed that a class of diagrams sums to zero on account 
of the single-field consistency relation [13, 33] (which is closely related to the Ward identity in 
(2.6); see [30] and Section 4), while others cancel because they are related by the non-linear 
transformation of £. 

These types of cancellations are reminiscent of those appearing in QED. For instance, con- 
sider photon-photon scattering. The leading-order diagram contains four external photons and 
a fermion loop connecting them. Each individual diagram, corresponding to a particular per- 
mutation of legs, is logarithmically divergent. However, the divergences exactly cancel when all 
diagrams are summed. In this case, the cancellation is, of course, a consequence of gauge invari- 
ance. To see this, consider the amplitude Mpuap, which by Lorentz invariance takes the following 
form 

M^vop = K(r]^ u r] ap + TJna-rjup + + finite terms . (2.14) 

A priori, the amplitude K could be divergent, but the Ward identity, p^Ai^ap = 0, forces it to 
be finite. This is an important result, since a divergence would have forced us to introduce a 
(A^A^) 2 counterterm, and consequently break gauge invariance [14]. In this paper, we will show 
that symmetries similarly protect correlation functions of curvature perturbations from getting 
a late-time evolution. 

Towards all orders. — At a qualitative level, it is easy to convince oneself that the dilatation 
symmetry implies constancy of £ to all orders in perturbation theory. We will ultimately agree 
with this intuition (see Section 3), but we would first like to point out where we feel that some 
details are missing. This may explain why some authors have not been convinced by these 
arguments. 

A general sentiment one encounters in the literature is that, because a constant £ mode can 
be removed by a gauge transformation, Cfc->o cannot be the source for a time-dependent solution. 
However, this argument appears somewhat circular since a time-dependent mode £(t) cannot be 
removed by such a transformation. On the other hand, one might have imagined that a time 
dependence of £ would require an operator equation of the form £ = ciQ + C2C 2 + • • • • The 
right-hand side of this equation is incompatible with the dilatation symmetry and is therefore 
forbidden to act as a source for £. However, why should such an operator equation be the only 
possibility? Moreover, a trivial counterexample to this logic is the case of a Goldstone boson, 
7r, which transforms as 7r \— > ir + 1. As in the case of £, a constant value of ir(x,t) = ttq is 
unphysical because it can be removed by a global transformation it \— > tt — ttq, which simply 
moves us between equivalent vacua. However, in flat space, the conclusion that 7Tfc_>o = as an 
operator is clearly false because quantum mechanics requires that [7^, 7r fc /] = i(27r) 3 5(k + k'). Of 
course, this counterexample isn't quite a fair analogy since in the case of £ we know that the 
modes become classical outside the horizon and freeze at tree level. Nevertheless, the example 
does illustrate that the argument has to involve more than symmetry alone. 

A more serious concern is that modes inside the horizon could induce a coherent effect on 
large scales that would cause a time dependence of the long- wavelength modes [6, 7]. These 
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short-scale modes are physical and cannot be removed by any symmetry. Hence, such coherent 
effects cannot be argued to vanish by symmetry alone. For example, time dependence could, in 
principle, arise from £(x,t) = cdi£d' l £(x,t), which is compatible with the dilatation symmetry. 
In momentum space, this becomes 

Ck = c J — 5» — CpCfc - p ' (2 - 15) 

which receives contributions from p > aH. The source term in eq. (2.15) is not suppressed as 
a — > oo since the momentum integral formally includes contributions from p — > oo. Of course, 
the treatment of these effects is complicated by the fact that the integral is UV divergent and 
needs to be regulated. One has to be careful that a bad choice of regulator doesn't introduce a 
spurious time dependence for (. 

Having described some of the subtleties involved in the conservation of £ at loop level, we 
will, in the next section, provide an all-orders symmetry-based proof for the time-independence 
of ("-correlators. 

3 A Non-Renormalization Theorem 

Our proof involves just a few relatively straightforward steps. First, we will prove that the 
modes always become classical outside the horizon (§3.1; see also [34]). We will show that this 
implies that any time evolution outside the horizon is described by an operator equation of the 
form 3 (§3.2) 

t(x,t) = Y,a n (t)C(x,t) + -- - . (3.1) 

n 

We will then use the dilatation symmetry to show that a n = (§3.3). Finally, we will show 
that the additional terms (• • • ) vanish at least like powers of a~ 2 and therefore can be ignored 
at late times. To understand this power law suppression requires a careful treatment of the 
renormalization of composite operators (see §3.4 and Appendix A). These terms include the 
effects of the modes inside the horizon that had been the concern of previous authors [6, 7]. 

Any no-go result is only as good as its assumptions. Let us therefore be clear about the 
assumptions that go into our proof: First, we will assume throughout that the theory is local and 
that the initial state is the Bunch-Davies vacuum. Second, we will only address loop corrections 
during inflation, such that the mode functions for the interaction pictures fields are roughly the 
de Sitter solutions. Finally, we will assume that any time-dependent couplings in the action 
for C scale at most like (loga(i)) r , for some finite r, and not as powers of a(t). The last two 
assumptions are mostly of technical nature and can probably be relaxed. However, even with 
these simplifying assumptions, our analysis is sufficiently general to cover the vast majority of 
inflationary models. We comment on ways to circumvent our theorem in §3.5. 

3 To avoid confusion, we will (in this section only) use a hat to denote quantum operators and reserve unhatted 
variables for c-numbers (such as the eigenvalues of (). 
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3.1 Locality and Classicality 

We begin by establishing the relation between locality of the theory and classicality of £ on 
superhorizon scales. We will define a mode £fc as being "classical" at late times, if it satisfies 

C = i\$lSL ^> , (3.2) 

<c 2 ><c 2 > 

where all the operators are evaluated at the same time. This definition of classicality implies that 
equal-time correlation functions of Cfe and/or can be rewritten in terms of classical stochastic 
variables, up to corrections that vanish as a — > oo — i.e. we can ignore all commutators at suf- 
ficiently late times. Moreover, eq. (3.2) assumes that the theory is approximately Gaussian, so 
that the power spectrum can be used to estimate of the size of any correlation function. 

Locality severely constrains the possible forms of equal-time commutators, like the one that 
appears in (3.2). In particular, the commutator of any pair of local operators 0\ and O2 must 
satisfy 

[6 1 (x,t),6 2 (y,t)] = , for x^y. (3.3) 
As a result, the commutator must be proportional (■>/— <?) _1 S(x — y) or derivatives therefore, 



[((x,t),((y,tj\ 



Y,Cn{t)d^\x,t)+[Y^d m {t)6 { [ n) {x,t))^d,+ 



S{X ~ y) , (3.4) 



where (n) and 6} m > are some basis of local scalar and vector operators, respectively. If the 
action for £ is time independent, then the coefficients in (3.4) must be time independent as well, 
i.e. c n (t) — > c n and d m [t) — > d m . Similarly, if the couplings in the action scale like (loga(t)) r 
for some finite r, then the coefficients in (3.4) are also logarithmic in a(t). In the limit a —> 00, 
we see that (3.4) therefore vanishes at least as a -3 , due to the overall factor of (V - 1 = a ~ 3 
required by diffeomorphism invar iance. 

To establish that the mode becomes classical in the sense of eq. (3.2), we now show that (£ 2 ) 
is bounded from below by a constant as a — > 00. First, let us insert a complete set of states into 
(2.11), 

(Qd(k) = ^(Qd|n)(n|C fe ) = (Qd|l)(l|Cfc) = \(^?5{k) , (3.5) 

n 

where we have rotated the basis of states such that (Qd\n) = S n i(Qd\l)- Similarly, we can insert 
the same set of states into the power spectrum of £ to find 

{CkC-k) = ^2\(Ck\n)\ 2 > |(Cfc|l)| 2 • (3.6) 

n 

Using Qd = and assuming 4 < 00, we must have |(C|1)| 2 > £ > where £ is a constant. 

We see that C — > as a — > 00 provided that vanishes more slowly than a" 3 . Recall that the goal 
of this section is to prove that £ vanishes at least as a -2 . Anything that violates our definition 
of classicality vanishes even faster. In that case, there is nothing for us to prove. 



4 This is essentially the assumption that a generalization of the Goldstone boson decay constant, /„., associated 
with Qd is finite. This is equivalent to demanding that ( is dynamical, i.e. has a finite kinetic term. 
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3.2 Operator Evolution 

We have proven that the modes of interest become classical at late times. If these were solutions 
to the classical equations of motion for a single degree of freedom, then they would be determined 
by two boundary conditions. For the problem at hand, one boundary condition is set by the choice 
of the Bunch-Davies vacuum and the other can be chosen to be the classical field configuration 
for £(cc, t) at a later time t. Therefore, given £(jc, t), the classical soultion for £(aj, t) is fixed. The 
purpose of this subsection is to make this statement precise, as an operator equation 

t(x,t) = f{C(x,t),dC(x, *),..•) , (3.7) 

for some functional /[(]■ 

Simultaneous eigenstates. — As in any quantum field theory describing a single degree of freedom, 
the operators £(x,t) form, at any time t, a complete set of commuting observables (one for each 
point in space). This has two important consequences [35]: 

• First, the eigenstates of these operators, \C(x,t)), are non-degenerate and form a complete 
basis of states on the Hilbert space. 5 

• Second, any operator 0(y,t) which commutes with £(x,t) is a function of £(x,t) alone, 
i.e. we have (D(y,t) = f[((y,t)], where / is a functional of £. 

Since at late times the commutator of £ and £ vanishes (see eq. (3.2)), we expect that, in the 
limit a(t) — > oo, the operator £ can be written as a function of £. Let us derive this result more 
formally. We start by defining the basis of eigenstates of £ as 

10 = |(Cl,C2, • • • ,&,"•)> , where ( a = ((x a ,t). (3.8) 

For clarity, we have used a discrete index to denote the spatial position. The state |£) is, by 
definition, an eigenstate of the operator £ a with eigenvalue ( a , i.e. 

Ca\C) = Ca|C) • (3.9) 
In this notation, the commutator (3.4) becomes 

[Lcb] = ■^j 5 «b+-' - ■ ( 3 - 10 ) 

where A a = ][] n c n (i)0i n ^(i) an d the ellipses denote terms that are suppressed by additional 
powers of a(t). Evaluating eq. (3.10) in the (-basis (3.8), we find 

(Cb - Cb)(C\L\C) = a tjf S ab + ■■■ ■ (3-11) 

5 The relation between these states and the classical solutions is most transparent in the Schrodinger picture, 
where we define a wavefunction for £(x, t), i.e. Vl/ff (as, t)]. For a single degree of freedom, the wavefunction satisfies 
a differential equation whose solution is determined by the initial state in the far past (e.g. Bunch-Davies) and the 
field configuration £(x, t) at late times t. See [13] for more details on the connection between in-in calculations 
and the Schrodinger representation. 
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This equation defines £ as an operator since it allows us to compute any matrix element by 
inserting a complete set of states. The r.h.s. of this equation scales at least as a -3 (up to log a 
corrections) and therefore the leading behavior is governed by the homogeneous solution, namely 

(C\L\C) « fa[CMC-C) + 0(a- 3 ) , (3.12) 

where S(C — C) — EL 'HCa — Ca) and 0(a~ 3 ) stands for operators whose correlation functions 
vanish as a — > oo (we make this more precise in Appendix A). 

Locality. — Next, let us see how locality constrains the form of the functional /o[C]- Recall that 
the conjugate momentum II satisfies the canonical commutation relation [^ a , lift] = i5 a b- Locality 
also requires that [^ a , flj,] oc Sab. Together with eq. (3.12), we then find 

<C| [L |C> = -^^p J (C " oc 5 ab . (3.13) 

As a result, f a [C] cannot depend explicitly on for b ^ a. Invariance under spatial translations 
furthermore implies that f a [(a] = /[Ca]- We have therefore established that 

(C\L\C)^f(Ca,dX a ,---)5(Ca-Ca). (3.14) 

Although this is a statement involving matrix elements in the C-basis, the result holds in any 
basis. To see this, note that (3.14) holds inside any correlation function: 

(L •••} = JdcJ dC(0|C)(Cl4lC)(CI •••) = (/[CJ •••) • (3.15) 

This proves that (3.14) is equivalent to the operator statement 

t(x,t) = f[((x,t)] , (3.16) 

where we have dropped the terms of order o -3 coming from the non-zero commutator (3.10). We 
will study the implications of this equation in the next two subsections. 

3.3 Constraints from Symmetry 

A basic property of any operator equation is that the two sides of the equation must transform 
in the same way under symmetries. In this subsection, we will show that the symmetries of C 
(see Section 2) severely constrain which operators are allow to appear on the r.h.s. of eq. (3.16). 

Perturbation expansion. — In perturbation theory, we usually consider situations where (C 2 ) "C 1. 
This corresponds to the requirement that the split of the metric into background and fluctuations 
is reliable. Given the small amplitude of fluctuations and the assumption of weak coupling, we 
can Taylor expand the r.h.s. of eq. (3.16) around C = 0, 

oo oo oo 

t(x, t) = J2 anWCix, t) + J2 Pm(t) {a- 2 e- 2 ^d 2 ax, t)) m + Y^ le(t) {<Fd£d£{x, t)) 1 + ■■■ . 

n=0 m=l l=\ 

(3.17) 



10 



Renormalization condition. — Since Q is a fluctuation, we require the expectation value of the 
l.h.s. of (3.17) to vanish, (£) =0. This fixes the coefficient of the unit operator, ao(t), in terms 
of the vacuum expectation values of the other local operators. Of course, we are always free to 
define (0) = as a renormalization condition for all local operators O ^ 1. In that case, the 
coefficient of the unit operator must vanish, cto(t) = 0. 

Dilatation symmetry. — Next, we consider the constraints imposed by the dilatation symmetry. 
Recall that i[Qd, C] = — 1 — x ' ®xCi which means that 

i[Qd,t\ =-x-d x t . (3.18) 

The higher-derivative operators in eq. (3.17) (i.e. those with coefficients /3 m , ye, etc.) have 
been arranged in such a way that i[Qd,0^] = —x ■ d x O^ . Therefore, any values of the 
coefficients f3 m and are consistent with the transformation of C( x i under Qd- The same is 
not true for the operators £ n , which transform as 

i [Qd, C] = -nC~ l ~ x ■ d x C ± -x ■ d x C ■ (3.19) 
We see that each individual term in the sum over Q n does not transform correctly to match the 
transformation of £. Furthermore, there is no way to choose the coefficients a n in such a way that 
the additional terms in the transformations of Q n cancel between terms. Therefore, consistency 
with the transformation under Qd requires that a n = for all n. 

Special conformal symmetry. — We can repeat the same analysis for the SCTs generated by Q l sc . 
From i[Q l sc , C] = — 2x* — 2x l (x ■ d x () + x 2 d l £, we infer that 

i[Qict\ =-2x l (x-d x t)+x 2 d l t . (3.20) 

Matching the transformation on the r.h.s. of (3.17) imposes non-trivial relations between the 
coefficients. For example, at second order in derivatives we have 

C(x,t) = a- 2 e- 2 ^{p 1 {t)d 2 C + li{tW 3 d l Cd ] C) + 0(d A ) . (3.21) 

Imposing that the transformations on both sides agree gives 

l{x,t) = p 1 (t)a- 2 e- 2 ^{d 2 C + \^d i Cd J C)+0{d i ) . (3.22) 

The special combination of operators on the r.h.s. should not be too surprising, since it is precisely 
the combination that appears in the three-dimensional Ricci scalar, 

U = -Aa- 2 e- 2i; (d 2 C + l&'diCdjC) . (3.23) 

Both the dilatation and the SCTs are continuously connected to a general, time-independent 
diffeomorphism on the spatial slice. Because £ transforms as a scalar under this group, the 
r.h.s. of (3.17) should be composed of invariants of the group. For this reason, our equation 
should take the form 

C{x,t) = j3 1 (t)K(x,t) + p 2 (t)n 2 (x,t) + fo^K^ix, + , (3.24) 

where IZij is the Ricci tensor on the spatial slice. The final step in our proof will be to show 
that all these terms vanish at least like powers of a~ 2 . Phrased in terms of curvatures, it seems 
intuitive that inflation should smooth out the spatial curvatures. Showing that this intuition 
survives quantum corrections will be the subject of the next subsection. 
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3.4 Renormalization of Composite Operators 

The right-hand-side of eq. (3.24) contains composite operators, i.e. products of fields evaluated 
at coincident points, which even in a Gaussian theory leads to divergences. One might worry 
that these divergences will affect the scaling behaviour of the operators at late times, i.e. change 
their a~ n suppression. In particular, the renormalization of composite operators is complicated 
by the tendency of operators to mix under renormalization [15, 16]. In order to complete the 
proof, we need to show that the scaling of the operators in (3.24) isn't drastically affected by 
renormalization. More precisely, we wish to show that if a composite operator 0(x,t) scales like 
a~ n in the free theory, any corrections in the interacting theory that scale like a~ m , where m is 
an integer with m < n, can be removed by a local redefinition of the operator, 

R (x,t) = 0(x,t) + 50(x,t) . (3.25) 

This allows us to define renormalized composite operators by local subtraction. By definition, 
these operators then all decay at least as powers of a~ 2 in correlation functions. 

An explicit demonstration of the renormalization of composite operators by local subtraction 
is rather technical. In this subsection, we therefore only show how the renormalization works in 
a concrete example (see also [36] for a related discussion). The dedicated reader can find the 
painful details for the most general cases in Appendix A. 

Example 

Consider a massless scalar field 6 in de Sitter space with interaction <p 3 . For purposes of illustration, 
we will present the renormalization of the composite operator 7 0(x) = (d 2 <fi / a 2 ) 2 (x) . 

Tree-level scaling. — Even in the free (or Gaussian) theory, this operator has a non-vanishing 
one-point function 



aA d 3 k k 4 ,, , 2 H 2 f aA k 3 dkf_ k 2 \ A 4 



(0)= 7^-M=jzi -zr- i + 



(2vr) 3 a 41 ^ 1 4tt 2 J a 4 \ ( aH ) 2 J 16tt 



H 2 + -A 2 
3 



(3.26) 



We have cut off the integral at fixed physical momentum A and used the Bunch-Davies mode 
function 

^ k (r) = ^=(l + ikr)e- ikT , (3.27) 

where r is conformal time. In accordance with our renormalization condition, we define a shifted 
operator with vanishing one-point function, 

Or = O - (O) . (3.28) 

Next, let us consider the two-point function of this operator (still in the Gaussian theory) 

(O R (x)O R (0)) = ^{{d 2 cj>{x)d 2 m)f > (3-29) 



6 Since the renormalization of composite operators is unrelated to the special symmetries satisfied by £, we have 
switched to a generic scalar field <f>. 

7 For notational simplicity, we will sometimes drop the time argument, i.e. O(x) means 0(x,t). 
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or 



j d 3 xe^ x (O R (x)O R (0)) = A | aA ^3g 4 |^| 2 |fc-Q| 4 |^- q | 2 



■q\k-q\ l + -r-FT^ l + 



2a 8 J (2vr) 3 *' ^' \ (aH) 2 J \ (aH)' 

1 H 4 k 5 ( 3 k 2 1 k A \ . 

1 + = t — H 7 — + (contact terms) . 



720tt 2 a 8 V 7(aH) 2 35 (ail) 4 y 

(3.30) 

We observe that the two-point function in the free theory scales as a -8 , as expected from the 
a -4 scaling of the operator. In the final line, we have dropped all terms that are analytic in 
k — e.g. (k 2 ) n , with n being a non-negative integer. If we Fourier transform such terms back to 
position space, they become contact terms — i.e. terms proportional to 8(x) — and therefore do 
not contribute to correlation functions at separated points. Notice that all terms proportional to 
the cutoff A are contact terms (as they should be for renormalized operators). 

One-loop correction. — Now consider the non-Gaussian correction to the cross-correlation 

1 f aA d 3 q 



where 



J d?xe ik - x {0 R {x)m) = ^f -^q 2 \k-q\ 2 {<t> q <j> k - q <i)- k )' , (3.31) 

(<Pq^h-q<t>-kY = i dfa(f) (H- mt (f)(t) q 4>k-q(f)-k(T)y + h.c. (3.32) 

J — oo 

The notation (•••)' denotes that an overall delta function has been omitted. Substituting the 
interaction Hamiltonian, 

Hint = 3~W J ' (3 ' 33) 

we get 

{</>q^k-q<f>-kY = ^<f>q(T)<Pi- q (T)<FL k (r)i j <j>' q (f) <t>' k - q (f) <f/_ k (f ) + h.c. , (3.34) 



where primes stand for derivatives with respect to conformal time. Since we are interested in the 
behavior as q — > oo, we keep only the leading terms in k, 

(MM . ig e - teig - M ijy^ + ,, , 

lff5ll 'l + 9 V + 2,VKl + 0f-)), (3-35) 



8 M 2 k q 5 

where K = q + \k — q\ + k. Note the importance of the Bunch-Davies vacuum in deriving 
eq. (3.35). In an excited state negative frequency modes would lead to contributions with K — > 
q — \ k — q\ + k ^ k. This would lead to extra inverse powers of k. Hence, we find 

1 1 H 5 1 f aA d 3 o 1 



/ 



1 1 1 / A 2 H 5 1A 4 H 3 1A 6 H\ „,,,„, 
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Notice that this cross-correlation scales as a~ 2 and not as a -4 (as we would naively expect from 
the scaling of the operator). The UV divergence has significantly affected the time dependence 
of the correlation function. This significant change in the scaling behavior of the operator would 
be a real problem if it weren't possible to remove the contribution by a local counterterm. 

Renormalization. — It is easy to see that the contribution in (3.36) can be removed by the following 
local operator 

1 / h?H 5 1A 4 # 3 1 A 6 H\ d 2 (t> 
d ° = ~^\~M 2 ~ + 2^d 2 ~ + ^^I 2 )^?' (3 ' 37) 



since 



|d 3 xe^(^)0(O)) = -fc 2 |^| 2 = -^ + o( (a ^) • (3.38) 

Moreover, we can also cancel higher powers of k in the expansion in (3.36). The first correction, of 
order k , is a pure contact term and therefore doesn't have to be removed explicitly. In fact, every 
even power, (k 2 ) n , where n is a non- negative integer, is a contact term and thus none of these 
terms contribute to correlation functions at separated points. This leaves the odd powers, k 2m ~ 3 , 
where m is a positive integer. It should be clear that all these terms can be removed by local 
operators of the form (d 2 ) m 4>. As a result, the contributions to the correlation function that lead 
to a physical scaling are associated with q ~ k <C aH (which we did not compute here). Clearly, 
all such contributions are suppressed by a -4 , as desired. In Appendix A, we will argue that 
defining renormalized operators by adding local counterterms is always possible in the Bunch- 
Davies vacuum. These operators have well-defined scaling behavior and are therefore suppressed 
at late times. In particular, operators with n derivatives vanish like (k/a) n . In Appendix A, we 
will also show that the renormalized operators satisfy the same symmetries as the bare operators. 
Higher-derivative composite operators therefore only make subleading contributions in eq. (3.24). 
This completes our proof. 



3.5 Summary of Results 

We have shown that 



' k 2 ' 

lim C k = + O ( -J ) . (3.39) 



Since this is an operator statement, it applies at all orders in the loop expansion. This means 
that any correlation function of will vanish as a — > oo. Equivalently, correlation functions of 
Cfe are time independent outside the horizon at all-loop order. 

We made four important assumptions in establishing this result: 

1. We assumed that we can transform to a gauge in which the scalar component of the metric, 
£, is the only propagating degree of freedom (in addition to gravitons). 

2. We assumed that the theory is local, in the sense that any pair of local operators 0\ and 
£?2 satisfies 

[6 1 (x,t),6 2 (y,t)] =0 , for x^y. (3.40) 

3. We assumed that couplings in the action for £ depends only logarithmically on the scale 
factor, i.e. X(t) oc (loga(i)) r , for some non-negative r. 
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4. We assumed the Bunch-Davies vacuum state. 
Using assumptions 1 - 3, we derived the following operator equation 

t(x,t) = f[C(x,t)] + 0(a- 3 ) . (3.41) 

Invariance under diffeomorphisms required that the lowest order terms in the derivative expansion 
are given by 

t(x,t) = p 1 (t)a- 2 e- 2 t(aPC + ±5 i Sd i Cd j C)+0(8 i ) . (3.42) 

This ensures that every operator in this series is suppressed by at least two derivatives. Finally, 
we showed that if assumption 4 holds, all operators containing derivatives are suppressed by 
factors of k/a and hence eq. (3.39) follows. 

Let us ask where our proof would fail if any of these assumptions were violated: 

1. In multi- field inflation, additional light scalars a are present in (-gauge. The proof that ( 
must satisfy an operator equation would still hold, but nothing would forbid terms of the 
form 

C(x,t) = K 1 {t)a(x,t) + K 2 (t)a 2 (x,t) + --- . (3.43) 

The fluctuations of a can be non-zero outside the horizon and therefore Cfc-s>o need not 
vanish. More dramatically, ref. [37] recently suggested an inflationary model (solid inflation) 
in which the adiabatic mode is completely absent and it isn't possible to go to the standard 
(-gauge. Our proof then doesn't apply. In fact, in solid inflation ( isn't conserved (even at 
tree level). 

2. Ref. [38] introduced an inflationary model (Khronon inflation) in which £ evolves as 

C fc (r) oc -^=e iakT , (3.44) 

where a is a ratio scales that will not matter here. We see that (fc_>o = i(k/a)<^k^07 which 
violates our eq. (3.39). However, one also finds that 

[Ur),Ck'(r)]=^5(k + k') , (3.45) 

which is non-local in real space. Khronon inflation therefore violates our locality assump- 
tion. 

3. In the model of ref. [39], the coefficient of the kinetic term Q 2 scales as a~ 6 (i), violating our 
assumption that couplings in the Lagrangian scale at most as loga(t). The authors of [39] 
then find solutions that scale as ( oc a 3 . This growing mode becomes classical and clearly 
satisfies an operator equation of the form ( ~ 3H(. Why is this equation not forbidden by 
the dilatation symmetry? First, we note that, due to the significant time dependence, the 
commutator scales as [Cfc(*)> Cfe'(*)] ~ ° 3 f° r a ~~ ^ 00 • ^ ms case > our operator equation 
takes the form ( ~ 3-ff( + O(aP). This is consistent with the dilatation symmetry because 
( i— > ( + A can be absorbed into 0(a ). 
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4. Assuming the Bunch-Davies vacuum was only important for the renormalization of opera- 
tors. In Appendix A, we show that large corrections to the scaling behavior of composite 
operators can be removed by a redefinition of the local operator. This renormalization 
procedure essentially requires that the only divergences in correlation functions arise from 
operators at coincident points. In some excited states, this is known not to be the case [40]. 
This is usually taken as a sign that the interacting theory is ill-defined. 

4 Relation to the Single-Field Consistency Relation 

In the previous section, we used both locality and symmetry to demonstrate that £ vanishes 
outside the horizon as a -2 . In the process, we understood the late-time scaling behavior of many 
other operators. In this section, we will see how this information is useful for understanding the 
behavior of correlation functions of Ck when there are large hierarchies between the momenta 
(i.e. for soft limits). Specifically, when the operator product expansion (OPE) applies, these 
correlation functions are determined in terms of the operators with the lowest scaling dimensions. 8 

4.1 Operator Product Expansion 

The OPE is a powerful tool for understanding quantum field theories in situations where the 
scaling behavior of operators is well understood [15]. The basic idea of the OPE is to replace 
a set of operators in the neighborhood of a point a; by a sum over local operators at x. When 
the operators can be organized according to their scaling dimensions (i.e. if one knows, for each 
operator, how many powers of the distance appear in correlation functions), then the leading 
contribution can be determined by the first few terms in the expansion. In the case of conformal 
field theories, one even understands the scaling behavior well enough to re-sum parts of the 
expansion. In the context of inflation, we would like to apply the OPE to correlation functions 
of £• On the surface, this doesn't look like a well-controlled procedure since ((x)((0) ~ log(|ac|) 
and therefore higher powers of £ are not suppressed in the OPE. However, in practice, the OPE 
is controlled by the smallness of (C 2 )- Furthermore, from the results of the previous sections (and 
Appendix A), we will be able to constrain the coefficient functions and/or the scaling behavior 
of each local operator. 

Consider the following OPE 

C(x)C(y) ^ £/ o (*_)0Oc + ) , (4.1) 
o 

where we defined 

x + = -(x + y) and X- = \x — y\ ■ (4-2) 
In Fourier space, this OPE reads 

C fc -|,U-| 9 ^ Efo(k)0- q . (4.3) 

o 

8 The utility of OPEs to describe soft limits of inflationary correlation functions has recently been emphasized 
by Kehagias and Riotto [41]. 
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The types of local operators O that should be included on the right-hand side depend on the field 
content of theory. Restricting to single-field inflation, the operators are composite operators made 
out of ( and its derivatives. The coefficient functions fo( u ) (° r their Fourier transforms fo(k)) are 
constrained by the symmetries of Section 2. Our arguments in Section 3 restrict the appearance 
of the operator £ in the OPE (in particular, up to corrections that vanish as a -3 , we can replace 
C using eq. (3.24)). Moreover, from the discussion in Appendix A, we know that higher-derivative 
composite operators are also suppressed by powers of a(t). The dominant operators in the OPE 
are therefore operators without derivatives 

COz)C(y) ^ £/„(*_)&(*+) + ■■■ • (4.4) 

n 

Acting n times with the dilatation charge, [Qd-, • • • ], on both sides of eq. (4.4), we find 

where £(x_) = (C(x)((y)) and P(k) = (CfeC-fc)'- The near-Gaussianity of the fluctuations, 
/nlC ^ 1) implies that it is often sufficient to keep only the n = 1 term in the sum 

|fc|»| g | dln(fc 3 P(fc)) 

C fe -| q C fe _| q ► / 1 (fc)C-, + -- = dlnfc C- q + - • (4.6) 

To capture subleading corrections, we should also consider derivative operators that have finite 
correlation functions in the limit a — > oo, such as 

C(aOC(y) ^ ■■■+gi(x-)xl[d l C](x + ) + 92(x^xl[d 2 ax + ) + ... , (4.7) 

where the functions g\(x—) and g%(x-) are at most logarithmic in x_. The Fourier transform of 
eq. (4.7) is 



SCI 



where the functions gi(k) and g^ik) scale as k~ . Acting with the charge of SCTs, [Q 
on both sides of eq. (4.8) implies that gi(k) = 0, which ensures that first subleading term is 
suppressed by q 2 /k 2 [12]. 

4.2 Single-Field Consistency Relation 

One is often interested in the coincident limit of correlation functions where two or more operators 
are brought close to each other (relative to the distances to other operator insertions). The OPE 
is a natural way to analyze this. For instance, let us consider the following limit of the three-point 
function lim x ^ y (((0)((x)((y)) . In Fourier space, this corresponds to the squeezed limit of the 
bispectrum, 

lim (CgCfc-IoC-fc-I-)' • (4-9) 

|q|«|fc| 2 Q 2 Q/ 

Maldacena showed that in single-field inflation this limit is fixed by the scale-dependence of the 
two-point function [13]. Here, we want to reproduce this result from the OPE. In fact, there is 
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not much left to do. We simply use the OPE (4.3) to replace the product of the high-momentum 
modes in (4.9), 

(C 9 C fe , ' ^ X>(*)fo0-,>' = +°(^))<^C-,> , + ... • (4.10) 

Terms that aren't shown explicitly in (4.10) are either suppressed by powers of a(t) or by the 
near-Gaussianity of the fluctuations. Hence, we find that the linear operators in the OPE (n = 1) 
lead to Maldacena's single-field consistency relation [13, 33] 

(C«Cfc-i,C- fc -i g ) ^(?) 

The vanishing of gi(fe) in (4.8) captures 9 the conformal consistency relation of Creminelli, Noreha 
and Simonovic [12], which ensures there is no 0{q ■ k) correction to (4.11). The subleading 
corrections (suppressed by q 2 /k 2 <C 1) come from operators like the last term in eq. (4.8). 

In [30], we related the single-field consistency relation to the Ward identity associated with 
the dilatation symmetry, eq. (2.6). Let us make a side remark addressed at readers familiar 
with that previous work. There we had to assume that 'multi-particle' states (i.e. states created 
by acting with several (^-operators on the vacuum) make a negligible contribution in single-field 
inflation. Here, we see that this assumption is equivalent to being able to truncate the OPE (4.3) 
at order n = 1. 

4.3 Violations of the Consistency Relation 

The conservation of £ and the consistency relation of the three-point function are closely related. 
In fact, both are consequences of the non-linearly realized dilatation symmetry. In §3.5, we list 
four critical assumptions on which our proof for the conservation of £ was based. It is interesting 
to see how violations of these assumptions map to proposed violations of the consistency relation: 

1. Of course, it is well-known that large squeezed limits are possible if £ isn't the only fluc- 
tuating degree of freedom (see e.g. [42, 43] for reviews of non-Gaussianity in multi-field 
inflation). In terms of the OPE, these additional fields may appear unsuppressed in (4.1). 
The coefficients of any such operator is not restricted by the non-linear symmetry of Q 
and therefore can be large. This is related, but not equivalent, to the violation of the 
conservation of £ we discussed in §3.5. 

2. In Khronon inflation [38] the consistency relation is still satisified, but subleading terms 
are now less suppressed (by q/k rather than q 2 /k 2 ). In fact, one could imagine variants of 
Khronon inflation that would even violate the consistency relation. It would be interesting 
to establish a more direct relation between non-locality and the scaling of the bispectrum 
in the squeezed limit. 

9 For higher n-point functions, the conformal consistency relation allows for terms that are linear in the soft 
external momenta [12]. This would arise from an OPE involving (n — 1) insertions of £ at separated points, but 
we will not consider such cases here. 



dln(k 6 P(k)) 



dink 



+ o 



k 2 



P(k) . 



(4.11) 
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3. Strongly time-dependent couplings (e.g. power law in a(t) rather than logarithmic) may 
violate the consistency relation directly through the non-conservation of £ [39] . In this case 
the operator £ scales as a n (for n > 0) outside the horizon and is therefore not suppressed 
in the OPE. This is also clear from the argument of Maldacena [13], which assumes that 
the k = mode can be removed by a rescaling of the coordinates. This is only true of 
time-independent solutions and therefore does not apply in this case. 

4. Some excited states have been found to violate 10 the consistency relation [44-47]. This 
can be understood as arising at times before the long mode has crossed the horizon, and 
hence derivative operators are not suppressed. However, since we directly apply the OPE 
at late times, it may not be clear where our argument breaks down. The resolution lies in 
the fact that the suppression of composite operators at late times (a — > oo) requires that 
we can remove any divergent result by a local counterterm. In these excited states, this is 
not the case [40]. However, this is an unphysical feature of these states that arises because 
they have infinite energy. One recovers the consistency relation in the k — > limit when 
restricting to finite energy states [44-47]. For an extensive discussion of these (and related) 
examples and their observational consequences, see [48]. 

5 Discussion 

In this paper, we proved that the super horizon conservation of the curvature perturbation ( in 
single-clock inflation holds as an operator statement. In the process, we developed techniques for 
understanding correlation functions of £ that did not require explicit use of perturbation theory. 
We used these insights to understand the operator product expansion of £ and its relation to the 
single-field consistency relation. 

There is reason to believe that the technical developments that we used to understand the 
conservation of £ may have applications to other problems. For example, we have not addressed 
the conservation of tensor modes. It should be clear that all constraints that followed from locality 
alone should apply equally to tensor perturbations. On the other hand, the transformation 
properties under large diffeomorphisms are quite different [11] and could lead to interesting results. 
Another application would be to eternal inflation. In this case, the fluctuations of £ are order one, 
which presents a challenge for using traditional perturbation theory techniques. However, our 
only result that made explicit use of perturbation theory was the scaling behavior of renormalized 
composite operators. For this reason, it is possible that some of our results will survive in the 
regime of eternal inflation. 

Finally, our primary concern was one type of infrared divergence of inflationary correlators, 
namely those that scale as loga(i). There are also infrared divergences that scale as logL, where 
L is a hard infrared cutoff on the comoving momenta. These types of divergences have been 
studied by many authors (see e.g. [49-56]) and it would be interesting to see if our understanding 
of the conservation of £ can shed any additional light on this other class of divergences. 



10 Technically speaking, the authors of [44-47] do not claim to violate the consistency relation in the limit k — > 0. 
However, if we removed the requirement that the states have finite energy, then a true violation can arise. Such 
states are essentially the same as those discussed in [40]. 
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Note added: When this paper was completed, ref. [57] appeared which also presents an all-orders 
proof for the conservation of £ on superhorizon scales. 
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A Renormalization of Composite Operators 

In quantum field theory in flat space, the Callan-Symanzik equation provides the connection 
between UV divergences and the scaling behavior of composite operators 0(x,t) [15, 16]. In 
perturbation theory, the scaling of composite operators is only corrected logarithmically, i.e. by 
small anomalous dimensions. This is manifest in renormalization schemes without power law 
divergences such as dimensional regularization. In schemes that allow for power law divergences, 
any power law corrections to the scaling can be removed by introducing local counterterms in 
the definition of renormalized operators Or. In this sense, only log corrections to scalings are 
physically meaningful. 

In this paper, we have been interested in the behavior of cosmological correlations functions 
in the limit a(t) — > oo. Here, the scale factor a(t) plays the role of an infrared regulator and the 
scaling with time t is controlled by the Hamiltonian and not the renormalization group. A priori, 
it is not obvious that there should be a relation between the results in flat space and in de Sitter 
space (although the two are mapped to each in the dS/CFT duality [58]). In this appendix, we 
will show explicitly that the intuition regarding anomalous dimensions in ordinary field theory 
will continue to hold for the scaling of cosmological correlation functions. In particular, we will 
show that the scaling behavior of composite operators is corrected in perturbation theory, at 
most, by loga(i). For example, suppose that, in the Gaussian theory, the two-point function of 
some composite operator 0(x,t) scales like a~ n , where n is some integer. We will show that 
at higher orders in perturbation theory, any contributions that scales like a~ m , where m is an 
integer with m < n, can be removed by a local redefinition of the operator, 



A.l Renormalizability by Local Counterterms 

In Section 3, we presented a specific example for the renormalization of a composite operator 
by local counterterms. This example was meant to be illustrative, but it does not establish that 
this procedure works at all orders in perturbation theory. In this section, we will put forward 
arguments to that effect. 

A. 1.1 Momentum Space Argument 

Perturbation theory is formulated most straightforwardly in momentum space, while locality is 
most manifest in position space. We will therefore present the argument twice, here in momentum 
and below in position space. In each case, we will first present the general strategy and then the 
details of the 'proof. 

Strategy. — Consider the composite operators 11 0(x,t) = Wi = i{oT 2 d 2 ) ri 4>{x : t). In momentum 
space, this becomes 



11 For notational simplicity, we will write most expressions for a specific composite operator, but our results hold 
for any local operator constructed from <j) and derivatives. 
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Regions of finite momentum in the integrand of (A. 2) scale like a~^ i<2ri in the limit a — > oo, 
relative to the correlation function of Y\2=i 0k ■■ Therefore, any contributions that do not scale 
like a~^ i2ri , must come from momentum configurations with ki — > oo and ki/a fixed (assuming 
that > 0). Because momentum is conserved, a least two of these momenta must diverge 

together. For purpose of illustration, let us consider the case where all the k^s diverge as a — > oo. 
We are interested in the scaling behavior of Ok inside of correlation functions, such as 

fim (O k <f> Pl ■ ■ ■ <f> p J = Em ( f[/ |^ E ^ " *K • • ■ K) • (A.3) 

i=l J v y i 



Our goal is to show that 

(ni- 2 

a,ki— too 



k' 

i im ( II Gi) 'V* i *(X) fei ~ fe )^Pi , "^m) = F (h,a)x]Ji> i 



i + o' p2 



(A.4) 



Let us explain why this is the desired result. First, notice that (A.4) factorizes into a function of 
the diverging momenta ki and a function of the finite momenta pj. This means that the integrals 
over the momenta ki in (A.3) will simply give a number C times the function of the momenta 
Pj. Second, the remaining function of pj is itself a correlation function of the fields (ft p . and some 
local operator at some lower order in perturbation theory. In the example above, the leading 
term would be the correlation function of f d 3 x e tk ' x (j) m (x) and <fi p ■ ■ ■ (j> p ■ On the other hand, 
if we were to find a factorized answer that contained higher inverse powers of pj (e.g. pj 5 ), then 
we would not be able to remove it by subtracting a local operator. Instead, we would need to 
subtract a non-local operator containing powers of d~ 2 . 

Our strategy will be to show that each Feynman diagram contributing to the above correla- 
tion function can be factorized into a sub-diagram containing the diverging momenta and one 
containing only finite momenta. After integrating over the momenta ki, we can simply replace 
the divergent sub-diagram with a local operator. 

In the following, we will generalize the standard arguments from flat space quantum field 
theory (see Weinberg, Vol. II, Ch. 20 [15]) to de Sitter space and the associated in-in correlation 
functions. Our argument will fall short of being a complete proof for technical reasons related to 
the regions of integration of the loop momenta. This is the same complication that arises in the 
standard arguments for validity of the OPE in flat space, like those in [15]. We will explain this 
in more detail at the end of this subsection. 

'Proof. — For concreteness, let us consider an (n + m)-point function of the form 

< K 2 ^) fcl • • • («~ 2 ^) fen ^ - O , (A.5) 

where all the operators are evaluated at some fixed time t. We are interested in the behavior as 
ki —> oo and a — > oo with k%ja fixed. One important feature of inflationary correlation functions is 
that the metric contains a trivial rescaling symmetry a — > \a and x — > \~ 1 x (and hence k — > Xk). 
Under this symmetry, local scalar operators transform as A and therefore their momentum space 
counterparts transform as A™ 3 . In general, this symmetry is broken in the action for £ because 
any explicit function of t can be rewritten in terms loga(t). In other words, we can't rescale a 
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and t independently after solving for the background. We will assume that the time dependence 
of the couplings is at most logarithmic in a(t) (power law in t, during inflation) such that the 
power law scaling, A~ 3 , is valid up to small corrections. Usually, this is required in order to 
preserve the near scale invariance of the fluctuations. As a result, we can think of the momentum 
dependence 12 of any field as being k~ 3 f(k/a)(l + \{n s — l)logfc + ••■), where k/a is counted 
as k°. 

Let us consider Feynman diagrams associated with the following in-in correlation function 

f dti ... dtrdH-^ih), . . . , [H int (t r ), {a~ 2 d 2 4>) ki ■ ■ ■ {a- 2 d 2 <f>) kn fa ■ ■■fa)))) . 

(A.6) 

We will focus on contributions from a diagram or sub-diagram, T, in which all internal momenta 
are of order k{ (for simplicity, we will take all the fcj's to be of comparable magnitude). Every 
external line contributes a factor 13 of k~ 3 or pj 3 . Since H- mi = J d 3 x a 3 (t)7i- m t(x, t) in position 
space, when written in momentum space, there is a momentum integral for every field and an 
overall momentum-conserving delta function. In terms of Feynman rules, this means that every 
internal line contributes a factor of k~ 3 for the contraction of the interaction picture fields and 
an integral J d 3 ki (recall that derivatives of the fields scale as k/a ~ k°). Since every vertex is 
associated with an insertion of H m t, each vertex introduces a momentum-conserving delta function 
and a time integral J dt'a 3 (t'), where we will count the factor of a 3 as k 3 f(k/a). In addition, 
there is a commutator associated with each vertex coming from the in-in expression (A.6). Each 
commutator with the fields 4> p is suppressed by p 3 /a 3 . These contributions can be ignored in 
the limit a — > oo. On the other hand, commutators acting on the internal lines do not affect the 
scaling since we are counting k/a as having scaling k°. 

Consider a general diagram with N vertices, I internal lines, external k lines and E p 
external p lines. The overall momentum scaling of the diagram is T ~ k D , where 

D = -3E k + 31 -3I + 3N - 3(N - 1) = -3E k + 3 . (A.7) 

Here, the term proportional to (N — 1) comes from extracting the overall momentum-conserving 
delta function. In writing this expression, we haven't been careful about the scaling with a(^) 
where tg (£ = 1, ■ ■ «,r) is the time appearing the l-th insertion of H[ nt . We are taking a(t) —> 
oo, but these integrals run over all values of ti, not just the far future. We have implicitly 
assumed that the i^-integrals receive their dominant contributions at late times, when a(ti) — > oo. 
Fortunately, this assumption is valid in the Bunch-Davies vacuum, because every internal lines 
comes with a factor of e - efc / a (*') from the ie prescription. As we take k —> oo, only contributions 
with a(ti) > ek aren't exponentially suppressed, and we may assume a(ti) — > oo inside F. 
Using these Feynman rules, we find that the diagram scales as 




(A.8) 

12 The mode functions in a scale-invariant theory are of the form <j>k{t) ~ k~ 3 ' 2 f{k/a(t)). Here, we are associating 
a fc~ 3//2 scaling with the creation and annihilation operators. 

13 For convenience, we will assume exact scale invariance of the interaction picture fields. 
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where K 3 is some product of three k^s and w,v are positive integers. The factors of \og(ki/H) 
can arise from the time dependence of the coupling constants. We do not include \og(p/H) 
scalings from time-dependent couplings, as they could only arise from early times. For Bunch- 
Davies initial conditions, these contributions are exponentially suppressed. Finally, we have 
included additional factors of loga(t) which, in principle, can arise from the time integration 
J*dt' ~ t ~ il" 1 log a [6]. 

We notice that eq. (A. 8) is precisely of the form of eq. (A. 4). We can therefore replace this 
whole diagram or sub-diagram by a local operator connected to the external lines with momenta 
pj times a function F(k{) (this is the operator production expansion in momentum space). When 
we define the correlation function in terms of the original local operator, we perform the integral 
over J d 3 k\ ■ ■ ■ d 3 k n 5(Yli &i ~~ YljPj)- Because the diagram factorizes, this is simply a number 
times a local operator. 

The higher-order terms in the p/k expansion come from Taylor expanding the internal lines 
in terms of pj. For example, an internal line might have momentum q = ki + p- ~ ki. Using 
rotation invariance, we can rewrite the dependence on q in terms of q ~ {kf + p 2 + 2fcj -pj) 1 ^ 2 . If 
we Taylor expand q in powers of pj , the only odd powers of pj appear in the combination ki ■ Pj . 
But odd powers of ki ■ pj vanish when we perform the angular integrals over the fcj's in (A. 2). 
With only even powers of pj surviving in the Taylor expansion, we can remove the entire series 
by adding derivatives inside the local operator. Therefore, all contributions to T that introduce 
power law changes in the scaling behavior of O can be removed by adding local counterterms. 

One may be concerned that the diagram V could be disconnected, in the sense that not all 
of the momenta ki are connected to each other via some path in V. Let us assume that this is 
true, i.e. let us assume that T is disconnected or connected only through soft internal lines. For 
momentum conservation to hold, this would imply that Y^i=i ~ ®(Pj)- This corresponds to 
a special momentum configuration where some subset of the diverging momenta separately sum 
to a finite momentum. For generic momenta this cannot arise and the diagram must therefore 
be connected. Since we will be integrating over ki, these non-generic points can be ignored. 

Loopholes and caveats. — The above argument seems very general, so it is worth highlighting 
situations where it could fail. First of all, the ie prescription of the Bunch-Davies vacuum 
was crucial for suppressing contributions at early times. For some excited states this may not 
be the case. This is consistent with the observation that some excited states in de Sitter are 
known not to be renormalizable by local counterterms [40]. In our language, these examples 
correspond to contributions when a(ti) ~ Pj/H, which would induce additional inverse powers 
of pj and therefore cannot be removed by local counter-terms. In extreme situations, one could 
also imagine compensating for the exponential suppression, e~ e «Jf, in the limit a — > 0, with 
exponential growth of the couplings of the form X(t) ~ exp(a _5 (i)) for 5 > 1. Having coupling 
grow this rapidly in the far past would give a large non-local contribution to late-time correlators. 
We have excluded such large time dependences throughout the paper. 

Finally, we want to stress that our argument falls short of a formal proof, as we have not 
been careful enough regarding the integration over loop momenta. Although we have treated all 
large momenta as order ki, there are always regions of integration where the internal momenta 
are much larger or much smaller than ki. One might worry that the result of performing and 
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regulating these integrals might somehow result in additional factors of ki /pj . The general scaling 
behavior of the internal lines makes it difficult to see how such contributions could arise, but we 
leave a complete investigation to future work. 

A. 1.2 Position Space Argument 

In translationally invariant theories, explicit calculations are often easier to perform in momentum 
space. However, locality is a fundamental property of most theories that is easier to understand 
in position space. For this reason, we will now explain how the renormalization of composite 
operators works in position space. The results will be less explicit than those of the previous 
section, but may be more intuitive. 

Strategy. — Consider the in-in master formula 

(0{x,t)0(y,t)} = (0\T e^-™^' Hint{t,) 0(x,t)0{y,t)T e-^-^^' Hintit,) |0> . (A.9) 
In perturbation theory this expression is evaluated by expanding in powers of 

J dt'H- mt (t') = J dt'd 3 x'a 3 (t')n hl t(x',t') (A.IO) 

and using contractions of the interaction picture fields. We then look at the regions of integration 
over the positions of Hi n t(,x', t') that are not suppressed by powers of a(t). We will show that 
these contributions arise from spacetime regions where some number of 'Hi n t(x' ,t') are within 
a sphere surrounding the composite operator 0(x,t) that is much smaller than the distances 
to other operators in the correlation function. By Taylor expanding Hmt(x' ,t') = Hmt(,x,t) + 

(x' — x) ■ d x Hmt(x, t') H , we can then treat 7imt(x' , t') as an operator at the point (x, t). To 

evaluate the behavior at coincident points, such as 0(x,t)T-Ci n t(x,t) (and derivatives therefore), 
we simply use the free field contractions of the fields that make up Hint and O. The result is 
therefore a new local operator with a divergent coefficient. The entire contribution can then be 
removed by adding a local counterterm to 0(x,t). 

The argument in position space is more subtle for the usual reason that perturbation theory 
is easier to implement in momentum space. For this reason, let us focus an a sightly simplified 
correlation function, namely the two-point function of the operator 0(x,t) = (a~ 2 d 2 4>) n (x , t) : 

/ dti ••• j 1 1 dt r {[H- m t{ti), ••• , [H- mt {t T ),0(x,t)0(y,t)}]}) . (A.ll) 

J —oo J — oo 

Despite the reduced complexity of the correlation function, it still is sufficient for our main goal. 
In particular, if the correlation function between (a~ 2 d 2 <fi) n (x , t) and any local operators is not 
suppressed (and cannot be removed by a local counterterm), then we can insert a complete set 
of states to find that ((a~ 2 d 2 4>) n (x, t)\n) is unsuppressed for some state |n). Inserting the same 
set of states in (A.ll) implies that there must be a similarly unsuppressed contribution to this 
two-point function. Each such contribution to (A.ll) takes the form \((a~ 2 d 2 cp) n (x,t)\n)\ 2 and 
therefore cannot be cancelled to make the final result vanish. Hence, it is sufficient to show that 
(A.ll) vanishes in the limit a — > oo to ensure that correlation functions of 0(x, t) with any local 
operator will be suppressed by powers of a -1 . 
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0(x,t) 0(y,t) 

. . t 




Figure 1: Composite operators 0(x,t) and 0(y,t) inserted at future infinity of de Sitter space. Locality 
requires that the operators are only influenced by sources that are inside their past light cones (shown in 
grey). The shaded regions near the operators indicate where an insertion of Hi nt (x' , t') yields an unsuppressed 
contribution to (0(x, t)0(y, t)). Divergences can therefore be removed by adding local counterterms. 

'Proof. — Let us first consider the contribution to (A. 11) from the insertion of 'Hmt(x r ,t r ) that 
includes a commutator acting on {d 2 <p) n at x or y. Locality requires that this commutator vanish 
when (x r ,t r ) is outside the past light-cone of (x,t) or (y,t). When t r ~ t, locality requires that 
the commutator is proportional to a delta function, 5(x r — x)0(x, t), or derivatives thereof (when 
acting on 0(x, £)). Any such contribution to the correlation function is manifestly local at x and 
can be removed by subtracting 0(x,t). 

Potentially dangerous contributions must come form points where t r <C t. In the Bunch- 
Davies vacuum, these contributions are localized on the past light cone of the operator 0(x,t), 
i.e. \x r — x\ ~ (a(t r )H) . This allows us to perform the integral J d 3 x r — > (a(t r )H)~ 2 f dfl. The 
factor of J dt r a 3 (t r ) from the measure in (A. 10) ensures that these contributions scale as a(t r ) 
in the limit a(t r ) — > 0. Naively, this suppression by a{t r ) in the measure could be compensated 
by the divergent contributions along the light-cone arising from the contractions of free fields. 
However, the ie prescription ensures that the only physical divergences come from operators at 
coincident points. Specifically, the ie prescription is equivalent to the analytic continuation of 
the correlators from the Euclidean sphere. These correlators are suppressed by powers of the 
Euclidean distance, which only vanishes at coincident points. The same will therefore be true 
of the Lorentzian correlations in the Bunch-Davies vacuum. Finally, the full expression must be 
invariant under the rescaling a — > Xa and x — > \~ l x. Due to the power law suppression along 
the light-cone 14 , our final result must be suppressed by (at least) a(t r )/a(t) for t r -C t. 

So far, we have only considered the contribution to the correlation function from a single 
insertion Wi n t(x r ,t r ). We found that, after integrating over x r and t r , only points near (x,t) or 
(y,t) contribute significantly (i.e. do not vanish as a — > oo). However, we should also consider 
what happens when the other interactions %\ n t{xi, ti) are included. From the above argument, the 
contribution from t r ~ t can still be removed by redefining the local operator, even in the presence 
of the additional interactions. This would suggest that the only contributions that cannot be 
removed are from points with t r <C t. However, we found that these points are suppressed by 
powers of a. The only way we would get a significant contribution would be if the integral over 
ti with insertion / H m t{ x ii'ti) would diverge as a positive power of a(t). However, for couplings 
that scale at most as loga(ij), the higher orders may diverge at most as (loga(i)) r [6]. This 

14 Otherwise, \x r — x\a(t r ) would be consistent with the rescaling symmetry and is unsuppressed. 
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completes our demonstration that composite operators can always be renormalized by adding 
local counterterms. 



A. 2 Symmetries of Renormalized Operators 

We conclude this appendix by showing that the basis of renormalized operators transforms under 
symmetries in the same way as the bare operators. 

First, let us review how to analyze the transformation of renormalized operators in the path 
integral formalism. Consider an action So that describes the full theory. To compute correlation 
functions, we deform the action by 

S = S + S n (J), (A. 12) 

where 

S n (J) = -J2j d 4 x^gJ n {x,t)( n {x,t) ■ (A.13) 
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We have added sources J n to the action such that correlation functions of £ n can be computed 
as 

"'>^o=^<"->L.' <A - 14) 

As far as the action is concerned, J n (x,t) is just like any other coupling. When we perturb 
in J n we will therefore find divergences that need to be removed by adding source-dependent 
counterterms 

Sr = S + S n (J) + 6 C S(J) . (A. 15) 

Because the theory is now finite, the correlation functions are also finite 



Sr,Jti—0 



The operators appearing in (A.16) are therefore the renormalized composite operators [C n ],R = Cr- 
In order to maintain the symmetries of the action, the sources J n can be given transformations 
under the symmetries such that S is invariant. Moreover, we can also choose these transformations 
to leave S n invariant. By taking the functional derivative of S n , we see that the operator jj- 
must transform in the same way as C n (x,t). Hence, as long as 5 C S(J) does not explicitly break 
the symmetries, then Sr is also invariant and therefore jj- = [C u ]r must transform in the same 
way as £ n . 

All of this applies equally to cosmological in-in correlation functions. In fact, in-in calcula- 
tions are just a special case of the above analysis, in which the (conformal) time integral goes 
from t = — oo(l — ie) to r = — oo(l + ie), while passing through r = 0. Because the symmetries 
of Section 2 are continuously connected to the group of diffeomorphisms, the counterterm action 
5 C S(J) would have to explicitly break diffeomorphism invariance to violate them. Using dimen- 
sional regularization and the results of the previous subsection, there is no need to use such a 
regulator. We therefore conclude that the transformation properties of renormalized operators 
follow from (2.4) and (2.5). 
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For the specific application in Section 3, there is, in fact, a more direct way to understand 
the symmetries of renormalized operators. In §3.3, we used symmetries to forbid the operators 
that transform non-linearly under the charge Qd- We then argued in §3.4 that the remaining 
operators vanish as a -2 after renormalization. The concern is that renormalization might mix 
these two groups of operators. For example, this would arise if On(k) = Ck — c(^k, where On 
scales as a -2 , C is some operator that transforms linearly under Qd and c is a constant. If this 
were the case, then we would find 



However, from eqs. (3.5) and (3.6), we see that any such mode has a constant contribution to its 
power spectrum. This violates the assumption that Oji(k) ~ a~ 2 , and therefore we must have 
c = 0. Repeating this argument for any other operators that vanish as a — > oo, we find that they 
all transform linearly under the dilatation symmetry. 




(A.17) 



28 



References 

[1] D. Wands, K. A. Malik, D. H. Lyth, and A. R. Liddle, "A New Approach to the Evolution of 
Cosmological Perturbations on Large Scales," Phys.Rev. D62 (2000) 043527, 
arXiv : astro-ph/0003278 [astro-ph] . 

[2] S. Weinberg, "Adiabatic Modes in Cosmology," Phys.Rev. D67 (2003) 123504, 
arXiv : astro-ph/0302326 [astro-ph] . 

[3] G. Rigopoulos and E. Shellard, "The Separate Universe Approach and the Evolution of Nonlinear 
Supcrhorizon Cosmological Perturbations," Phys.Rev. D68 (2003) 123518, 
arXiv : astro-ph/0306620 [astro-ph] . 

[4] D. H. Lyth, K. A. Malik, and M. Sasaki, "A General Proof of the Conservation of the Curvature 
Perturbation," JCAP 0505 (2005) 004, arXiv: astro-ph/0411220 [astro-ph]. 

[5] D. Langlois and F. Vernizzi, "Conserved Non-Linear Quantities in Cosmology," Phys.Rev. D72 
(2005) 103501, arXiv:astro-ph/0509078 [astro-ph]. 

[6] S. Weinberg, "Quantum Contributions to Cosmological Correlations. II. Can These Corrections 
Become Large?," Phys.Rev. D74 (2006) 023508, arXiv:hep-th/0605244 [hep-th] . 

[7] S. Weinberg, "Quantum Contributions to Cosmological Correlations," Phys.Rev. D72 (2005) 
043514, arXiv:hep-th/0506236 [hep-th]. 

[8] L. Senatore and M. Zaldarriaga, "On Loops in Inflation," JHEP 1012 (2010) 008, 
arXiv: 0912. 2734 [hep-th]. 

[9] P. Creminelli, M. A. Luty, A. Nicolis, and L. Senatore, "Starting the Universe: Stable Violation of 
the Null Energy Condition and Non-Standard Cosmologies," JHEP 0612 (2006) 080, 
arXiv : hep-th/0606090 [hep-th] . 

[10] C. Cheung, P. Creminelli, A. L. Fitzpatrick, J. Kaplan, and L. Senatore, "The Effective Field 
Theory of Inflation," JHEP 0803 (2008) 014, arXiv: 0709 . 0293 [hep-th]. 

[11] K. Hinterbichler, L. Hui, and J. Khoury, "Conformal Symmetries of Adiabatic Modes in 
Cosmology," JCAP 1208 (2012) 017, arXiv: 1203.6351 [hep-th]. 

[12] P. Creminelli, J. Norena, and M. Simonovic, "Conformal Consistency Relations for Single-Field 
Inflation," JCAP 1207 (2012) 052, arXiv: 1203.4595 [hep-th]. 

[13] J. M. Maldacena, "Non-Gaussian Features of Primordial Fluctuations in Single Field Inflationary 
Models," JHEP 0305 (2003) 013, arXiv:astro-ph/0210603 [astro-ph]. 

[14] S. Weinberg, The Quantum Theory of Fields, vol. 1. Cambridge University Press, Cambridge, 2005. 

[15] S. Weinberg, The Quantum Theory of Fields, vol. 2. Cambridge University Press, Cambridge, 2005. 

[16] J. Collins, Renormalization : An Introduction to Renormalization, The Renormalization Group, and 
The Operator- Product Expansion. Cambridge University Press, Cambridge, 1984. 

[17] L. Senatore and M. Zaldarriaga, "The Effective Field Theory of Multifield Inflation," JHEP 1204 
(2012) 024, arXiv: 1009.2093 [hep-th]. 

[18] L. Senatore and M. Zaldarriaga, "A Naturally Large Four-Point Function in Single Field Inflation," 
JCAP 1101 (2011) 003, arXiv : 1004 . 1201 [hep-th] . 

[19] D. Baumann and D. Green, "Equilateral Non-Gaussianity and New Physics on the Horizon," JCAP 
1109 (2011) 014, arXiv: 1102. 5343 [hep-th]. 



29 



[20] D. Baumann and D. Green, "Signatures of Supersymmetry from the Early Universe," Phys.Rev. 
D85 (2012) 103520, arXiv: 1109.0292 [hep-th]. 

[21] D. Baumann and D. Green, "A Field Range Bound for General Single-Field Inflation," J CAP 1205 
(2012) 017, arXiv: 1111.3040 [hep-th]. 

[22] S. R. Behbahani and D. Green, "Collective Symmetry Breaking and Resonant Non-Gaussianity," 
arXiv: 1207.2779 [hep-th]. 

[23] A. Nicolis and F. Piazza, "Spontaneous Symmetry Probing," JEEP 1206 (2012) 025, 
arXiv: 1112.5174 [hep-th]. 

[24] S. R. Behbahani, A. Dymarsky, M. Mirbabayi, and L. Senatore, "(Small) Resonant 

Non-Gaussianities: Signatures of a Discrete Shift Symmetry in the Effective Field Theory of 
Inflation," arXiv : 1111 . 3373 [hep-th]. 

[25] R. Gwyn, G. A. Palma, M. Sakellariadou, and S. Sypsas, "Effective Field Theory of Weakly 
Coupled Inflationary Models," arXiv : 1210 . 3020 [hep-th] . 

[26] D. Lopez Nacir, R. A. Porto, L. Senatore, and M. Zaldarriaga, "Dissipative Effects in the Effective 
Field Theory of Inflation," JEEP 1201 (2012) 075, arXiv: 1109.4192 [hep-th]. 

[27] D. Lopez Nacir, R. A. Porto, and M. Zaldarriaga, "The Consistency Condition for the Three-Point 
Function in Dissipative Single-Clock Inflation," JCAP 1209 (2012) 004, arXiv : 1206 . 7083 
[hep-th] . 

[28] J. M. Bardeen, P. J. Steinhardt, and M. S. Turner, "Spontaneous Creation of Almost Scale-Free 
Density Perturbations in an Inflationary Universe," Phys.Rev. D28 (1983) 679. 

[29] D. Salopek and J. Bond, "Nonlinear Evolution of Long Wavelength Metric Fluctuations in 
Inflationary Models," Phys.Rev. D42 (1990) 3936-3962. 

[30] V. Assassi, D. Baumann, and D. Green, "On Soft Limits of Inflationary Correlation Functions," 
arXiv: 1204.4207 [hep-th]. 

[31] E. Kahya, V. Oncmli, and R. Woodard, "The Zeta-Zeta Correlator Is Time Dependent," Phys.Lett. 
B694 (2010) 101-107, arXiv: 1006.3999 [astro-ph.CO] . 

[32] G. L. Pimentel, L. Senatore, and M. Zaldarriaga, "On Loops in Inflation III: Time Independence of 
Zeta in Single Clock Inflation," JEEP 1207 (2012) 166, arXiv: 1203.6651 [hep-th] . 

[33] P. Creminelli and M. Zaldarriaga, "Single-Field Consistency Relation for the Three-Point 
Function," JCAP 0410 (2004) 006, arXiv: astro-ph/0407059 [astro-ph] . 

[34] D. H. Lyth and D. Seery, "Classicality of the Primordial Perturbations," Phys.Lett. B662 (2008) 
309 313, arXiv:astro-ph/0607647 [astro-ph]. 

[35] P. Dirac, The Principles of Quantum Mechanics. Clarendon Press, 1958. 

[36] D. Harlow and D. Stanford, "Operator Dictionaries and Wave Functions in AdS/CFT and 
dS/CFT," arXiv: 1104. 2621 [hep-th]. 

[37] S. Endlich, A. Nicolis, and J. Wang, "Solid Inflation," arXiv : 1210 . 0569 [hep-th]. 

[38] P. Creminelli, J. Norena, M. Pena, and M. Simonovic, "Khronon Inflation," arXiv : 1206 . 1083 
[hep-th] . 

[39] M. H. Namjoo, H. Firouzjahi, and M. Sasaki, "Violation of Non-Gaussianity Consistency Relation 
in a Single-Field Inflationary Model," arXiv: 1210.3692 [astro-ph. CD] . 

[40] T. Banks and L. Mannelli, "De Sitter Vacua, Renormalization and Locality," Phys.Rev. D67 (2003) 



30 



065009, arXiv:hep-th/0209113 [hep-th] . 

[41] A. Kehagias and A. Riotto, "Operator Product Expansion of Inflationary Correlators and 

Conformal Symmetry of de Sitter," Nucl.Phys. B864 (2012) 492-529, arXiv : 1205 . 1523 [hep-th]. 

[42] D. Wands, "Multiple-Field Inflation," Lect.Notes Phys. 738 (2008) 275-304, 
arXiv : astro-ph/0702187 [astro-ph] . 

[43] D. Langlois, "Non-Gaussianities from Isocurvature Modes," arXiv: 1209 . 0330 [astro-ph. CO] . 

[44] I. Agullo and L. Parker, "Non-Gaussianities and the Stimulated Creation of Quanta in the 
Inflationary Universe," Phys.Rev. D83 (2011) 063526, arXiv : 1010 . 5766 [astro-ph. CO] . 

[45] I. Agullo and L. Parker, "Stimulated Creation of Quanta during Inflation and the Observable 
Universe," Gen.Rel.Grav. 43 (2011) 2541-2545, arXiv: 1106.4240 [astro-ph. CO] . 

[46] J. Ganc, "Calculating the Local- Type Jnl for Slow-Roll Inflation with a Non- Vacuum Initial 
State," Phys.Rev. D84 (2011) 063514, arXiv: 1104. 0244 [astro-ph. CO] . 

[47] D. Chialva, "Signatures of Very High Energy Physics in the Squeezed Limit of the Bispectrum 

(Violation of Maldacena's Condition)," JCAP 1210 (2012) 037, arXiv : 1108 .4203 [astro-ph. CO] . 

[48] R. Flauger, D. Green, and R. Porto, "to appear," arXiv: 1211 .xxxx [hep-th]. 

[49] D. Seery, "Infrared Effects in Inflationary Correlation Functions," Class. Quant. Grav. 27 (2010) 
124005, arXiv : 1005 . 1649 [astro-ph . CO] . 

[50] S. B. Giddings and M. S. Sloth, "Semiclassical Relations and IR Effects in de Sitter and Slow-Roll 
Spacetimes," JCAP 1101 (2011) 023, arXiv: 1005. 1056 [hep-th]. 

[51] C. T. Byrnes, M. Gerstenlauer, A. Hebecker, S. Nurmi, and G. Tasinato, "Inflationary Infrared 
Divergences: Geometry of the Reheating Surface versus SN Formalism," JCAP 1008 (2010) 006, 
arXiv: 1005.3307 [hep-th]. 

[52] M. Gerstenlauer, A. Hebecker, and G. Tasinato, "Inflationary Correlation Functions without 
Infrared Divergences," JCAP 1106 (2011) 021, arXiv: 1102.0560 [astro-ph. CO] . 

[53] S. B. Giddings and M. S. Sloth, "Cosmological Observables, IR Growth of Fluctuations, and 
Scale-Dependent Anisotropics," Phys.Rev. D84 (2011) 063528, arXiv: 1104. 0002 [hep-th]. 

[54] S. B. Giddings and M. S. Sloth, "Fluctuating Geometries, q-Obscrvables, and Infrared Growth in 
Inflationary Spacetimes," arXiv: 1109. 1000 [hep-th]. 

[55] S. Miao and R. Woodard, "Issues Concerning Loop Corrections to the Primordial Power Spectra," 
JCAP 1207 (2012) 008, arXiv : 1204 . 1784 [astro-ph . CO] . 

[56] L. Senatore and M. Zaldarriaga, "On Loops in Inflation II: IR Effects in Single-Clock Inflation," 
arXiv: 1203.6354 [hep-th]. 

[57] L. Senatore and M. Zaldarriaga, "The Constancy of Zeta in Single-Clock Inflation At All Loops," 
arXiv: 1210.6048 [hep-th]. 

[58] A. Strominger, "The dS/CFT Correspondence," J HEP QUO (2001) 034, arXiv:hep-th/0106113 
[hep-th] . 



31 



